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Abstract 

Chaotic dynamics of the dripping faucet was investigated both exper- 
imentaUy and theoretically. We measured continuous change in drop po- 
sition and velocity using a high-speed camera. Continuous trajectories of 
a low-dimensional chaotic attractor were reconstructed from these data, 
which was not previously obtained but predicted in our fluid dynamic 
simulation. From the simulation, we further obtained an approximate po- 
tential function with only two variables, the drop mass and its position 
of the center of mass. The potential landscape helps one to understand 
intuitively how the dripping dynamics can exhibit low-dimensional chaos. 

It is weW known that the beat of a dripping water faucet is not always 
regular and exhibits complex behaviour including chaos [Q. Chaotic dripping 
was originally suggested by Rossler |Q and experimentally confirmed by Shaw 
and his collaborators ||^ . Since then, many experimental and theoretical studies 
have established the dripping faucet as a sort of paradigm for chaotic systems 

i i I I B 111 6 111' B @ & Most previous studies 
have involved measuring the time interval T„ between successive drips, because 
the dripping time is easily measured using a drop-counter apparatus ^, ^, ^ 
H H' 0- "^^^ time intervals are then plotted in pairs (T„, T„+i) for each n to 
give a return map. Because the return maps typically appear low dimensional (~ 
one-dimensional), the behavior is often described by a simple dynamical model 
composed of a variable mass and a spring ^, 0, ^ |4[ |l^ . In this mass- 

spring model, a mass point, whose mass increases linearly with time at a given 
flow rate Q, oscillates with a fixed value of the spring constant fc; and a part Am 
of the total mass m is removed when the spring extension exceeds a threshold, 
which describes the breakup of a drop. Although the model exhibits chaotic 
return maps similar to those obtained experimentally, its empirical nature means 
that it does not provide a unified explanation for the complex behavior of the 
real dripping faucets, typically seen in bifurcation diagrams (plotting of T„ vs. 
the control parameter Q). Thus, the connection between the low-dimensional 
dynamical system and a presumably infinite-dimensional fiuid dynamical system 
remains elusive. 

Recently, we carried out fiuid dynamic computations (FDC) based on a new 
algorithm involving Lagrangian description and succeeded in reproducing not 
only the time-dependent shapes of the drops, but also various characteristics of 
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nonlinear dynamical systems, such as period-doubling, intermittency, hystere- 
sis, observed in experiments of dripping faucets p6[ |. A detailed analysis of our 
FDC made it possible to improve the mass-spring model by using more realistic 
approximations including taking into account the mass dependence of the spring 
constant, and the necking process Q. This improved mass-spring model, de- 
scribed only by three variables: mass m, its position zq and the velocity zq, 
closely described the dynamical behavior exhibited by the FDC. In particular, 
qualitatively good agreement among bifurcation diagrams (especially the repeat- 
ing structure which reflects the oscillation frequency of each drop), obtained by 
(i) the experiment jl^, (ii) the FDC jl^ and (iii) the improved mass-spring 
model , over a wide range of the control parameter Q provided a theoretical 
model of the basic low-dimcnsional structure inherent in the system, which is 
presented here. 

An experimental chaotic attractor is shown in Fig. |^. In addition to the 
discrete-time variable T„, the mass m(t) of the liquid suspended by the nozzle, 
the position of the center of mass zg(^) and velocity zq were also measured. 
As seen in Fig. |](b), ZQ{t) oscillates because the surface tension works as a 
restoring force. Similar oscillation was observed in both our FDC and improved 
mass-spring model Note that the frequency of the oscillation during each 
time interval T„ is determined by the flow rate. At Q = 0.24 g/s, zq oscillates 
6 ± 1/2 times before the breaking up and the phase of the oscillation in the 
plane {zg^zq) fluctuates chaotically at breakup. A return map of {zpeak.n} for 
every fourth peak of zait) is presented in Fig. |l|(c). The return map of T„ (Fig. 
0(a)) looks like an entangled string and thus the map function T„+i = G{Tn) 
cannot be single- valued. In contrast, the map function 2;peak,n+i — Fizpciik,n) is 
approximately single-valued, although broadened by 'noise'. Low-dimensional 
multi- valued map functions of T„ have often been obtained in both experiments 
@' ^^"^ numerical simulations |16, For the improved mass-spring 



model ]17| , the map function mr,n+i = H{m-c^n) (wr,n being the remnant mass 
just after the n-th breakup), as well as the map function F(zpeak,rt) is single- 
valued, while G(Tn) is generally multi-valued. This is because the mapping 
of the time interval T„ is not topologically conjugate to that of the dynamical 
variable mr^„. The FDC results also showed that Him-c^n) is approximately 
single-valued even when G(Tn) is multi-valued. Thus, T„ does not directly 
specify the state of the system despite its easy measurement. 

Measurement of the continuous-time variables {zoit), ZGit),m{t)} made it 
possible to visualize the trajectory of the chaotic attractor in a continuous state 
space for the first time. The projection of the attractor in the plane {zq,zq) 
corresponding to the chaotic motion in Fig. ^ is shown in Fig. |(top). On av- 
erage, zcit) oscfllated six times during each time interval r„ (see Figs. |l|(b) 
and ||(top)). The cross section of the attractor in Fig. |^(a) looks nearly one- 
dimensional and explicitly shows that the drop formation is well characterized 
by a few state variables. Figure ^(b) is a blow-up of the region S in which tra- 
jectories starting with slightly different remnant masses separate rapidly from 
each other. The transition from the oscillating process to the 'necking' process 
(i.e., stretching of the attractor) occurs in this region S |l|]. In Fig. |(bottom) 
the FDC results are presented including profiles of the evolving drops (Fig. 
^(d)), which show qualitatively good agreement with the experimental observa- 
tions. The cross section of the theoretical attractor (Fig. ^(c)) also exhibits the 
low-dimensionality. An attractor in the space {zQ,ZQ,m) with a similar spiral 



2 



structure has already been obtained using the improved mass-spring model ||T^ . 
These results strongly support the use of the three state variables: zq, zq and 
m to describe the drop motion. 

The motion of the drop is subjected to gravitational force and surface ten- 
sion. Since the surface energy depends on the shape of the drop, the total 
potential energy U (i.e., gravitational plus surface) should be a function of the 
many degrees of freedom of the liquid. Our FDC show, however, that approxi- 
mating U as a function of two variables, m and zq, yields a conceptually clear 
picture for the basic low-dimensional structure of the system. Cross sections of 
U/m, the potential energy per unit mass, at several fixed values of m for the 
chaotic attractor shown in Fig. ||(bottom) are presented in Fig. ||(a). Poincare 
cross sections at the same values of m as Fig. ||(a) are given in Fig. ||(b). Al- 
though the cross sections in Fig. ^(a) look two-folded, if they are approximated 
as single- valued functions, then a sheet of the potential surface U{m, ZQ)/m is 
obtained as shown in Fig. ^(c). The surface is characterized by a U-shaped 
valley and a ridge which converge as m increases and have totally merged when 
TO = TOcrit, the maximum mass of the static stable shape. We have also carried 
out a numerical experiment based on FDC with Q = by fixing the value of 
TO. Various artificial shapes of the drops were used as initial conditions, and it 
was found that the point (to, zq) rapidly converged onto a well-defined potential 
surface U{m, ZG)/m similar to that shown in Fig. ||(c) |21 



The potential surface for the closed system (i.e., the pendant drop for (5 = 0) 
is useful for understanding the dynamics of the open system {Q > 0). The 
minimum and a maximum of U(m, ZQ)/m for the Q = system correspond 
to a stable fixed point (static equilibrium shape) and a saddle point (unstable 
equilibrium shape, see Fig. ||(c)), respectively. These are indicated by red points 
in the top two panels of Fig. ||(b). When to — TOcrit {— 0.144 g), the stable fixed 
point disappears, and the necking process begins leading to the pinching off of 
the drop. When to < merit, the trajectories of the drops oscillate in the valley of 
U{m, ZQ)/m, which corresponds to the spiral structure of the attractor in Fig. 
2. The flow of the Q = system as depicted in Fig. ^(b), also suggests that the 
state point of the drop causes the spiral orbits in the phase space (zq, ZG,m). 
The chaotic dynamics, especially the stretching and folding of the attractor, 
can be interpreted as follows. The stable fixed point approaches the saddle 
point as to increases. The state point thus visits the vicinity of the saddle point 
after the spiral motion, which causes the instability of the orbit, or equivalently, 
the stretching of the attractor (as in Fig. ||(b)). The Poincare section of the 
chaotic attractor (Fig. ||(b)), which rotates around the stable fixed point with 
increasing to, is gradually crushed and asymptotes to a single line through the 
necking process. 

The above analysis leads to the motion of the dripping faucet being described 
in terms of the approximate potential function U as 

dU{m,ZG)/m (7 + Q) . . ^ 

ZG^ Z ,TO = Q, (1) 

OZQ TO 

where 7 is the damping constant. This equation describes the motion of a 
particle with unit mass in 2-dimensional space (to, zq), where the particle is 
forced to have a constant velocity m = Q in the direction of the TO-axis. Our 
improved mass-spring model corresponds to a local quadratic approximation for 
the potential near the bottom zq = zq of the valley in which the TO-dependence 
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of the spring constant fc(m) = d'^U{'m, ZQ)/dzQ has been included jl^. The 
present description is thus a natural extension of the improved mass-spring 
model. As we have seen, however, the picture of the global landscape of the 
potential is much clearer than the improved mass-spring model for intuitively 
understanding the dynamics of the system. If the flow rate Q is small enough so 
that the initial oscillation after the breakup is damped, the particle (equivalently 
the state point of the drop) goes along the bottom of the valley as m increases. 
At TO = TOcrit, where the valley merges into the ridge, the particle's motion loses 
its stability (start of the necking process) and rapidly approaches the breakup 
point (on the broken red line in Fig. ||(c)). Since the instability always begins 
at the same point m = TOcrit in this case, the size of each drop is uniform. If, 
on the other hand, the oscillation of the drop affects the start of necking, the 
dripping motion exhibits a variety of flow rate-dependent periodic and chaotic 
patterns. Note that the particle generally gets over the ridge before to reaches 
TOcrit when it is oscillating. At a flow rate resulting in a chaotic motion, two 
trajectories starting with slightly different to values (different remnant masses) 
are initially close to each other. As to increases, however, one trajectory may 
get over the ridge at a certain m value, while the other remains in the valley 
owing to a small difference in zq (two red lines in Fig. ^(c)). 

Similar potential landscapes are expected for other subjects involving drop 
formation like ink-jet printing [ p2[ and atomization | |2^ . The potential surface 
describing the motion in the attractor was confirmed by FDC to have essentially 
the same topology even for high viscosity liquids and for microfluidic systems 
(drop size ~ pi). Hence, the present picture of a particle moving along a low- 
dimensional potential surface is relevant to modeling drop formation in general. 

We are grateful to Professor W. S. Price and Professor A. L. Salas-Brito for 
careful reading of the manuscript. 
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Figure 1: Experiments on the dripping nozzle (7 mm inner diameter and 10 mm 
outer diameter). The flow rate Q = 0.24 g/s (~ 2 drips/s on average), (a) Re- 
turn map of dripping time intervals, where the time intervals in pairs (T„, T^+i) 
for each n are plotted. The time series {T„} was measured using a drop-counter 
apparatus wherein a laser beam directed at a detector is interrupted whenever 
a drop crosses the beam, (b) Oscillation of the position zq of the center of mass 
of the drop under the nozzle, zq was estimated from the shape of the drop 
using the digitised image recorded by a high-speed video camera every 1 /500 s. 
(c) Return map of the fourth peak ^;peak,n of the oscillating zq, where Zpeak,n is 
indicated by an arrow in (b). 
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Figure 2: The top panel is the experimentaUy observed strange attractor under 
the same experimental conditions as in Fig. where the position zq and the 
velocity zq of the drop under the nozzle are plotted at every 1/500 s. (a) 
Poincare cross section of the attractor (plotting of zq vs. zq) at t = 0.3 s 
measured from each breakup moment; (b) Enlarged view of the region S in 
which the attractor undergoes stretching and folding. The bottom panel is 
numerical simulations, where water drops falling from a nozzle (7 mm diameter) 
are simulated at the flow rate Q = 0.32 g/s. We have used the same parameter 
values as in Ref. (c) Poincare section of the attractor {zq vs. zq) at 

m = 0.13 g; (d) Drop deformation. 
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Figure 3: (a) Fluid dynamic simulation (7 mm diameter, the flow rate Q — 0.32 
g/s). Cross section of U /m for the attractor is plotted vs. zq on a constant mass 
surface, where U /m is the total potential (= surface tension plus gravitation) 
energy per unit mass as a function of m (drop mass in units of gram in each 
panel) and zq (position of the center of mass). The critical mass merit = 0.144 
g. (b) Black dots represent the Poincare section of the attractor, plotting of zq 
vs. 2Gj on the same constant mass surface as in (a) (obtained from the same 
simulation). The red solid curves indicate the flow in the state space, which were 
drawn by referring to the results obtained from another simulation with fixed 
mass (i.e., Q — Q). The red points in the top two panels are fixed points for the 
Q = system, (c) Schematic view of a typical potential surface, obtained by 
approximating U /m in (a) by a single valued function of zq and to. The drop 
shapes (light blue) in various state points (zg,to) on the potential surface are 
indicated. The red lines represent the trajectories starting with slightly different 
masses; one oscillates in the potential valley, while the other gets over the ridge 
because of a slight difference of the velocity zq. The breakup occurs on the red 
dashed line. 
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